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ABSTRACT 



Context. The iron lines at 630.15 and 630.25 nm are often used to determine the physical conditions of the solar photosphere. A 

common approach is to invert them simultaneously under the Milne-Eddington approximation. The same thermodynamic parameters 

are employed for the two lines, except for their opacities, which are assumed to have a constant ratio. 

Aims. We aim at investigating the validity of this assumption, since the two lines are not exactly the same. 

Methods. We use magnetohydrodynamic simulations of the quiet Sun to examine the behavior of the ME thermodynamic parameters 

and their influence on the retrieval of vector magnetic fields and flow velocities. 

Results. Our analysis shows that the two lines can be coupled and inverted simultaneously using the same thermodynamic parameters 

and a constant opacity ratio. The inversion of two lines is significantly more accurate than single-line inversions because of the larger 

number of observables. 
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1. Introduction 

. Both solar physics and radiative transfer owe many things to 

■ the Milne-Eddington (ME) approximation. Being an incredibly 
' simplistic description of spectral line formation, it provides use- 
„ ful hints to understand the behavior of spectral lines as formed 

■ in solar and stellar atmospheres. It also offers a key diagnos- 
' tics to infer the physical conditions of the plasma. This is par- 
ticularly true when magnetic fields are present. The first solu- 
tion of the radiative transfer equation in the presence of a mag- 
netic field was derived adopting the ME assumption that all 

■ physic al quantities relevant to line forni ation are constant with 
I depth (IUnnolll956 l : iRachkovskvl 1 1 962allbb . Under these condi- 
. tions, the solution is analytic and, therefore, by simply vary- 
ing the model parameters, one gets a handle on the behavior of 
Stokes line profiles. Similarly, the analytic character of the ME 

. solution allows perturbative analyses li ke those performed by 
iLandi Degl'Innocenti & Landolfil (1 19831). who studied the influ- 
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[ ence o f velocity gradients, or lOrozco Suarez & Del Toro Iniesta I 
(l2007l) . who calculated the sensitivities of ME Stokes pro- 
" files to the various model parameters. Inversion codes of 
' the radiative transfer equation are diagnostic tools that be- 
come simpler under the M E hypo thesis. Since the p ioneer - 



ing work by iHarvev et aP (Il972b and lAueretal.] (| 1977|) 



and a fter i mprovements b y iLandolfi & Landi Degl'Innocentil 
(Il982h and iLandolfi et all (ll984|), a number of ME inver- 
sion c odes have been develope d . These include th e HAO 
code (ISkunianich & LitesI [19871 iLites & Skumanict] 
MELAN IE (ISocas-Navarroll2001h. HEL IX (iLagg et all 
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MILOS ((Oro zco Sua rez & Del Toro Iniesta 1 12007|), and VFISV 
(iBorrero et al.. .20 10). 

Strictly speaking, the ME model is applicable to just one 
spectral line. The reason is that the so-called thermodynamic 
parameters of the model are meant to characterize the behavior 



of the specific line under consideration. The line-to-continuum 
opacity ratio, rjo, the Doppler width of the line, A/Id, and the 
damping parameter, a, govern the shape of the Stokes profiles 
(i.e., they are the parameters of the Voigt and Faraday- Voigt 
functions). In turn, the source-fuction terms, 5*0 and 5"!, con- 
trol the continuum level, the line depression, and the Stokes am- 
plitudes. However, many investigations are based upon the si- 
multaneous ME inversion of spectral line pairs, like the well- 
known Fei doublet at 630 nm (e.g., Lites et al. 1993 ) or the 
Mgi b lines at 517.2 and 518.3 nm (iLites et al.iri988h . These 
inversions are reported to use no extra free parameters as com- 
pared to single-line inversions, based on the similarities between 
the two lines belonging to the same multiplet: the opacities for 
each line are specified in the ratio of their respective oscillator 
strengths, while the Doppler widths and the damping parameters 
are assumed to be identical for the two lines. 5*0 and S i are also 
the same for both lines since they lie very close in wavelength. 

This strategy has been widely used for more than 20 years 
to analyze the Fe 1 630 nm measurements ta ken with instruments 



such as the Advanced Stokes Polarimeter (JElmore et al.lll992]) 



or the spectropolarimeter aboard the Hinode s atellite (JLites et al.l 
120071: iKosugiet al.l2007l:lTsuneta et aT]l2008h . The simultaneous 
inversion has been s hown to provide better results than single- 
line inversions (e.g.. iLites et al.l[T994l) . The better performance 
is easy to understand: if the observables reproduced by the inver- 
sion are doubled, the results can be expected to be more accurate 
(at least by a factor V2). However, the two lines do not have ex- 
actly the same atomic parameters (see Table[T]). This makes them 
to be formed at slightly diflfe rent heights, where diflTerent ph ysi- 
cal conditions may exit (e.g. lMartmez Gonzalez et al1l2006h . In 
view of these diflTerences and of the large excursions of 770, a, 
and A/1d in the real solar photosphere, both horizontally and ver- 
tically, one may wonder whether the simultaneous ME inversion 
of the two lines is valid and what the limitations are. General 
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Table 1. Atomic data for the Fe i 630 nm lines. 



Ao (nm) 


;r(eV) 


log^/ 


Transition 


alai 


(T 


^eff 


630.1501 
630.2494 


3.654 
3.686 


-0.75 
-1.236 


5P2 - 5Z)2 
5Pi - 5Do 


0.243 
0.240 


840.5 
856.8 


1.67 

2.5 



Notes. Shown are the central wavelength of the transition, Xq, the exci- 
tation potential of the lower atomic level, ;^, the multiplicity of the lower 
level times the oscillator strength, log g/, the collisional broadening pa- 
rameters, a and cr (in units of Bohr's radius ao), and the effective Lande 
factor of the line, geff- The logg/ values have been derived from a fit 
to the solar spectrum using a two-component model of the quiet Sun 
(iBorrero & Bellot Rubioll2002h . 



inversion codes not relying on the ME approximation perform 
exact line transfer calculations, so they are able to invert the 
two lines without inconsistencies. The purpose of the present 
Research Note is to check the ME case using the excellent test 
bench offered by modern magnetohydrodynamic (MHD) simu- 
lations. 



2. Influence of the ME thermodynamic parameters 

Trade-offs among the ME thermodynamic parameters have been 
reported and found not to be very important for the infer- 
ence of vec tor magnetic fields and l ine-of-sight (LOS) ve- 
lociti es (e.g.. iLites & Skumanich|[l99Ql: IWestendorp Plaza et al.l 
19981). An explanation of this phe n omen on has been given by 
Orozco Suarez & Del Toro Iniesta I (l2007l) . Therefore, the first 
question we should answer is whether or not the thermodynamic 
parameters of the lines need to be the same to obtain similar LOS 
velocities, vlos, magnetic field strengths, B, inclinations, y, and 
azimuths, ^, from their individual analysis. 

To investigate this issue, we synthesize realistic Stokes pro- 
files for the two Fei lines at 630 nm using simulations per- 
formed with the MPS/University of Chicago Radiative MHD 
code. This code solves the MHD equations for compressible 
and partially ionized plasma s. Fu rther information abo ut the 
simulations can be found in IVog ler (2003) and Vogle retaD 
(l2005l) . The snapshot used here belongs to a mixed-polarity sim- 
ulation run with an average magnetic field strength of {B) ^ 
140 G at log T500 = -1. We generate the Stokes spectra 
of the two lines from the MHD m odels using the SIR code 
(iRuiz Cobo & Del Toro Iniestalll992h . as explained by Orozco 
Suarez et al. (2010). We then invert the two lines separately with 
the MILOS code. The inversion is carried out assuming a one- 
component model atmosphere (magnetic filling factor unity) and 
zero macroturbulent velocity. No noise is added to the Stokes 
profiles. 

Figure[T]shows scatter plots of A/Id, ci, Sq, and 770 as inferred 
from the inversion of the two lines. The results for S 1 are very 
similar to those for 5*0, with less scatter. In this and other fig- 
ures, the gray colors inform about the pixel density, with black 
meaning larger values. Over-plotted are dashed lines represent- 
ing one-to-one correspondences, except for 7/0 where the ratio 
^0,2 / ^0,1 = 0.327 is indicated (1 and 2 stand for the 630.15 and 
the 630.25 lines, respectively; see below). The plots include all 
the pixels in the simulation snapshot (288 x 288), independently 
of the polarization signal. Obviously, the thermodynamic param- 
eters are far from being the same for the two lines. The scatter 
is large for A/Id and 5*0, although they show a linear correlation 
with a slope close to unity. In the case of 7/0 and a the scatter is 
dramatic; differences of up to a factor 3 for a and of more than 
one order of magnitude for 7/0 can be seen. As a matter of fact. 
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Fig. 1. A/Id, ^, So, and rjo parameters from the inversion of the 
Fe I line at 630.25 nm vs those obtained from the 630. 15 nm line. 
The dashed lines represent one-to-one correspondences. In the 
bottom right panel the dashed line represents 7/0,2/^0,1 = 0.327. 
Both lines are inverted separately. 
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Fig. 2. LOS velocity, magnetic field strength, inclination, and az- 
imuth from the inversion of the 630.25 nm line vs those from the 
inversion of the line at 630.15 nm. Both lines are inverted sepa- 
rately. 



the 7/0 values obtained for the 630.15 nm line span the full range 
of variation of the line-to-contin uum opacity ratio in real sola r 
atmospheres (see, e.g.. Fig. 11 of lWestendorp Plaza et al.|[T998l) . 

Do these discrepancies between the thermodynamic param- 
eters affect the magnetic and velocity inferences? The answer to 
this question is negative as can be seen in Fig. [2l where vlos, ^, 
7, and if are displayed as inferred from the individual inversion 
of the 630.15 and 630.25 nm lines. Despite the strong scatter in 
the thermodynamics, the results from both lines are remarkably 
similar. Therefore, we have to conclude that although the ther- 
modynamic parameters of ME inversions may have little mean- 
ing, it is possible to establish approximate relations between the 
parameters of the two lines for use in simultaneous inversions, 
as we shall see in the next section. 
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3. Understanding the simultaneous inversion 

Since iron is mostly ionized and local thermodynamic equilib- 
rium conditions prevail in the solar photosphere, the dependence 
of the logarithmic derivative of t]o on temperature (the most im- 
portant quantity in line formation) is linear with the excitation 
potential;^ and the ionization potential /, and quadratic with the 
temperature. A/Id is proportional to the square root of tempera- 
ture and does not depend on the atomic parameters of the line. 
Assuming van der Waals broadening for the calculation of the 
damping coefficient, an intricate (but weak) dependence of a on 
X, along with a proportionalit y with the 7 /lQ-th power of tem- 
perature, is obtained (see, e.g. jGrayll2008l) . 

Neglecting the differences induced by the slightly different 
excitation potentials of both lines, the a and A/Id ratios are pro- 
portional to /lo,2/^o,i- Thus, the damping and Doppler width can 
safely be assumed equal for both lines, as is also the case for 5*0 
and S 1 because the Planck function should not present signifi- 
cant differences in such a short wavelength interval, t/o is pro- 
portional to the square of the central wavelength of the line and 
to its gf value. Therefore, ^0,2/^0,1 ~ (^2 f2)/(gi /i), where g is 
the multiplicity of the lower level and / the oscillator strength. 
Note that this ratio is independent of the temperature or other 
thermodynamic quantities, so it should not change with height 
in the atmosphere even if the opacity varies by orders of magni- 
tude. Using the parameters of Table 1, the opacity ratio for the 
Fei 630 nm lines turns out to be 7/0,2 / ^0,1 = 0.327. This is the 
value implemented in MILOS. The HAO code and MELANIE 
use similar ratio^. 

To check the validity of this estimate, we invert the 
630.25 nm line again but now forcing 5*0, Si, A/Id, and a to 
be equal to those obtained from the previous inversion of the 
630.15 nm line. The remaining model parameters are allowed to 
vary freely. Figure [3] shows the 7/0 values retrieved from the in- 
version. The dashed line represents the ratio 770,2 / ^0,1 = 0.327. 
The solid line corresponds to a multi-polynomial fity = a[l]x^ -\- 
a[2]x^ -F a[3]x with a = [lx lO'^, -9.3 x 10"^ 0.5] for 770,1 < 38 
and y = b[l]x-\- b[2] with b = [0.23, 4. 1] for 770,1 > 38. Note that 
the theoretical ratio provides a fair description of the relation- 
ship between the two 770 values in the range where most of the 
points are located (to stress the differences, the inset zooms in on 
the boxed area). Since the exact values of the thermodynamical 
parameters are not very important for the determination of the 
magnetic field vector and the LOS velocity, we conclude that it 
is safe to use a constant opacity ratio to invert the two lines si- 
multaneously without increasing the number of free parameters. 

The final consistency proof is shown in Fig. ID where the 
parameters obtained from the inversion of Fei 630.25 nm are 
plotted against those coming from the simultaneous inversion 
of the two lines as coupled through their theoretical 770 ratio. 
The scatter is very low for B and vlos and somewhat larger for 
y and (p, but still much smaller that that of Fig. O This sug- 
gests that the accuracy of analyses based on one single line (e.g., 
iBommier et al.l 120091) could be improved by adding the other 
line. 
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vestigated in this paper. Our work has been supported by the Spanish MICINN 



^ The opacity ratio depends basically upon /, which is known with 
limited precision and varies from one source to the next. This causes 
an uncertainty in the th eoretical opacity ratio. For example, the labo- 
ratory measurements of lBardetal. II 19911 give log gf = -0.718 for Fei 
630.15 nm, rather than the -0.75 specified in Table 1. With this oscil- 
lator strength, the theoretical ratio would be 0.301 (closer to the values 
retrieved from the inversion, but only at the high 770 end). 




Fig. 3. 770 values retrieved from the inversion of the Fe I lines at 
630.15 and 630.25 nm. The dashed line represents the theoret- 
ical ratio 770,2/^0,1 = 0.327. The solid line stands for a multi- 
polynomial fit. The inset zooms in the dotted box. 
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Fig. 4. IMagnetic field strength, inclination, azimuth, and LOS 
velocity from the the inversion of the 630.25 nm line vs those 
from the simultaneous inversion of the two lines. The latter in- 
version is done assuming the same model atmosphere (total of 9 
free parameters) for the two lines, with their 770 values coupled. 
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